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1.   Introduction. 

The  underlying  problem  to  be  treated  here  is  that  of  inter- 
polation of  scattered  data.   Given  irregularly  spaced  points 
(xk'vk'^k)»  k  =  lf«fN/  with  the  assumption  that  the  (xk,yk) 
values  do  not  all  lie  on  a  straight  line,  the  problem  is  to 
construct  a  smooth  function,  F(x,y),  which  takes  on  the  values  fk 
at  the  points  (*k,yk).   This  problem  has  received  considerable 
attention  over  the  last  few  years,  and  a  great  deal  has  been 
learned.   Several  surveys  of  the  field  are  available,  including 
Schumaker  (1976),  Barnhill  (1977),  and  Franke  (1982).   The  first 
two  treat  a  wider  class  of  approximations.   The  latter  contains 
the  results  of  testing  a  large  number  of  algorithms  for  the 
problem.   Despite  the  existence  of  many  methods  for  construction 
of  interpolating  surfaces,  a  single  scheme  suitable  for  a  wide 
range  of  dispositions  of  the  independent  variable  data  and  for 
functions  which  have  large  gradients  implied  by  the  data  has  not 
yet  been  devised. 

The  problem  of  overshoot  by  surfaces  in  the  presence  of 
large  gradients  has  been  previously  addressed  by  Nielson  and 
Franke  (1984)  from  the  point  of  view  of  generalizing  the 
univariate  spline  under  tension.   There  the  construction  was 
based  on  an  underlying  triangulation  of  the  independent  variable 
data,  and  the  surface  reflected  this  as  tension  was  increased. 
While  potentially  useful  for  user  control  of  the  surface,  such  a 
scheme  clearly  does  not  model  a  physical  process,  as  does  the 
univariate  spline  under  tension  (see  Schweikert  (1966)  and 
references  in  Nielson  and  Franke).   Another  appoach  to  modeling 
such  surfaces,  based  on  the  assumption  that  the  surface  actually 


has  a  discontinuity  has  been  investigated  by  Franke  and  Nielsori 
(1983).  Automatic  detection  of  probable  discontinui tes  was  not 
considered. 

The  present  paper  is  an  attempt  to  investigate  a  scheme  for 
representation  of  surfaces  which  incorporates  the  concept  of 
tension  in  a  way  which  models  a  physical  process.   The  deflection 
of  a  thin  plate  under  the  imposition  of  point  loads  in  the  lat- 
eral direction  and  mid-plane  forces  (tension)  will  be  developed. 
The  theory  of  interpolation  of  scattered  data  by  thin  plate 
splines  (TPS)  under  point  loads  is  well  known.   Harder  and 
Desmarais  (1972)  developed  the  equation  from  engineering  consi- 
derations.  This  was  followed  by  an  elegant  mathematical  analysis 
and  generalization  by  Duchon  (1975,1976,1977)  and  Meinguet 
(1979,1979a).   Further  generalization  to  smoothing  scattered  data 
was  given  by  Wahba  and  Wendelberger  (1980). 

The  construction  of  the  interpolation  function,  once  the 
appropriate  basis  functions  are  determined,  will  parallel  that  of 
TPS  and  other  methods  based  on  the  association  of  a  basis  func- 
tion with  each  data  point.   Let  the  function  Bk(x,y)  be 
associated  with  the  point  Uk,yk).   The  inclusion  of  some  poly- 
nomial terms  in  the  approximation  may  arise  either  in  a  natural 
way  from  the  mathematics  or  physics  of  the  problem,  or  from  a 
desire  to  incorporate  polynomial  precision  into  the  method.   With 
linear  terms  the  approximation  then  becomes 

N 
F(x,y)  =  Y  AkBk(x,y)  +  a  +  bx  +  cy   . 

The   coefficients    in    the   approximation  are  determined   by   requiring 


that  the  function  interpolate,  and  that  certain  conditions  on  the 
coefficients  be  satisfied.   The  latter  conditions  are  related  to 
the  polynomial  terms,  and  imply  the  method  has  some  polynomial 
precision  (for  linear  functions,  above).   For  TPS,  linear  func- 
tions are  in  the  kernel  of  the  functional  being  minimized.   In 
addition,  the  constraint  equations  for  TPS  may  be  thought  of  as 
"equi libr ium"  conditions  on  the  plate;   the  sums  of  the  forces 
and  moments  must  be  zero.   In  the  end,  the  coefficients  are 
determined  by  the  linear  system 


AkBk(xi,yi)  +  a  +  bxi  +  cyi  =  f^      ,  i  =  l,...,N 


5 

k  =  l 
N 

I 

N 

1 

N 

y  A*yk  = 
k^l 


Ak  ■  0 


Akxk  =  0 


(1) 


■  0 


The  development  of  the  equation  of  deflection  of  a  thin 
plate  under  the  influence  of  point  loads  and  mid-plane  stresses 
will  be  developed  from  the  engineering  point  of  view  in  the  next 
section.   A  discussion  of  the  properties  of  the  basis  functions, 
and  their  behavior  as  tension  goes  to  zero  or  becomes  large  will 
be  given  in  Section  3.   Some  examples  showing  the  effectiveness 
of  the  scheme  in  controlling  overshoot  and  discussion  of  some 
problems  of  the  scheme  will  be  given  in  the  Section  4.   Finally, 
some  speculations  concerning  a  more  elegant  mathemat ization  of 
the  method  will  be  discussed. 


2.   Basis  Functions  for  Thin  Plate  Splines  with  Tension 

The  development  of  the  basis  functions  for  thin  plate 
splines  with  tension  (TPST)  will  approximately  parallel  that  of 
Harder  and  Desmarais  for  surface  splines  (as  they  called  thin 
plate  splines).   We  begin  with  the  equation  of  a  thin  plate 
subjected  to  lateral  loads  and  mid-plane  forces.   This  equation 
is  found  in  books  on  plates  and  shells,  for  example,  Timoshenko 
and  Woinowsky-Krieger  (1959, p.  379).   In  its  general  form  the 
equation  is 

VW    _(q  +  Nx_2-   +  Ny_      Nxyg^T   )    ' 

where  W  is  the  lateral  deflection,  q  is  the  lateral  load,  Nx,  Ny, 
and  N    are  the  mid-plane  forces,  and  D  depends  on  the  properties 
of  the  plate  material.  Setting  Ny  =  Nx  and  Nxv  =  0,  one  obtains 
an  equation  of  the  form 

V4W  -  «2V2W  =  p   , 

where    is  a  tension  parameter.   The  solution  of  the  equation 
will  be  determined  by  a  two  step  process.   Let  V  =  V2w,  then  we 
have   V2V  -  «  2V  =  p. 

The  deflection  of  the  plate,  Wa,  under  a  point  load  at  the 
origin  (the  Green's  function)  will  now  be  developed.   Radial 
symmetry  is  assumed  and  converting  the  equation  to  polar  co- 
ordinates under  this  assumption  yields 

1  d_(  dV}  _  a2v  .  6 

r  dr  v  dr  ' 

where  6  is  the  unit  impulse  functional.   This  is  a  modified 


Bessel  equation  of  order  zero,  and  the  two  solutions  of  the 
homogeneous  equation  are  IQ(ar)  and  KQ(ar).   Among  the  boundary 
conditions  to  be  applied  are  that  the  second  derivatives  of  Wa 
tend  to  zero  as  r — »°°   .   Thus  IQ  is  rejected  and  the  solution  is 
found  to  be  -  (2ir)   KQ  (ar ) .   The  factor  -(2-rr)"1  is  obtained  to 
satisfy  the  jump  condition  at  the  origin.   Details  are  given  in 
Stakgold  (1979, p.  77).   The  remaining  part  of  the  problem  now  is 
to  solve  the  equation 


1  d  ,  dwrt, 

—  ^— ( r a )  = 

r    dr   dr 


-(2Tr)-1K0(ar) 


While  no  explicit  elementary  function  is  obtained,  two  integra- 
tions with  application  of  the  condition  of  finiteness  at  the 

origin  and  an  arbitrary  value,  C,  at  the  origin  yield 

r    t 

Wa(r)    =   -(2-n)'1  f  t'1  JsKQ(as)ds6t   +  C  .  (2) 

o  o 

The  lack  of  an  explicit  representation  of  the  function  is  not  a 
hindrance  either  theoretically  or  computationally.   The  details 
of  computing  with  this  basis  function  will  be  discussed  in  the 
next  section. 

By  the  process  of  superposition  of  point  loads,  and  inclu- 
sion of  an  "equilibrium  condition"  for  the  plate,  one  obtains  the 
representation 

N 


F(x,y)  =  ^AkWa(rk)  +  a 


(3) 


k  =  l 


where  r^  =  (x-xk)2  +  (y-yk)2.   The  coefficients  are  obtained  by 
solution  of  the  system  of  equations  corresponding  to  Eq.  (1),  but 
without  the  x  and  y  terms  in  the  approximation,  and  hence  without 


the  last  two  conditions.   The  constraint  equation  requires  that 
the  sum  of  the  loads  (coefficients,  Ak),  be  equal  to  zero.   As 
will  be  seen  in  the  next  section,  this  condition  could  be  derived 
by  requiring  that  the  surface  become  "flat"  as  r  —  »  .   Because  of 
the  inclusion  of  the  constant  term  in  the  approximation,  one  can 
take  the  value  of  C  in  (2)  to  be  anything  convenient.   In  my 
computations  I  took  C  =  0,  but  another  value  could  lead  to  better 
numerical  stability  in  the  solution  of  the  system  for  the  Ak  and 
a,  especially  for  large  values  of  tension.   This  was  not  investigate 


3.   Properties  of  the  Approximation 

We  first  note  that  the  lack  of  an  elementary  representation 
of  the  basis  function  is  not  a  serious  problem.   The  function  can 
be  easily  approximated,  either  by  numerical  approximation  of  the 
integral  formulation,  above,  or  by  considering  the  function  as 
one  of  the  components  of  the  solution  of  a  pair  of  ordinary 
differential  equations.   In  the  numerical  examples  given  in  a 
later  section,  the  latter  approach  was  used.   Using  the  ODE 
solver,  DVERK  from  the  IMSL  library,  a  table  of  values  of  the 
function  was  obtained  over  a  suitable  range  for  the  problem  at 
hand,  and  with  a  spacing  which  allowed  linear  interpolation  for 
values  at  intermediate  points  to  the  desired  accuracy  (in  my 
case,  to  essentially  single  precision  accuracy  on  an  IBM  compu- 
ter).  To  illustrate  the  behavior  of  the  function  for  various 
values  of  the  tension  parameter.  Figure  1  gives  plots  of  the 
function  for  several  parameter  values.   Here,  each  function  has 
been  normalized  to  have  value  1  at  r  -   2.   Multiplication  of  a 


basis  function  by  a  constant  has  no  effect  on  the  overall  approx- 
imation since  the  constant  is  accounted  for  in  the  coefficient. 
As  the  tension  parameter  is  increased,  the  basis  function  Wa 
becomes  more  and  more  "pointed"  at  the  origin.   Because  the 
equation  of  the  thin  plate  with  tension  approaches  the  membrane 
equation  (Laplace's  equation)  as  tension  gets  large,  and  because 
the  Green's  function  for  the  membrane  equation  is   log  r,  this  is 
not  surprising. 

According  to  the  differential  equation  the  overall  approxi- 
mation (if  not  the  individual  basis  functions)  should  approach 
the  TPS  when  the  tension  is  decreased  to  zero.   However,  this 
cannot  be  the  case  since  the  TPS  includes  linear  terms  in  the 
approximation,  while  the  TPST  includes  only  a  constant  term. 
Nevertheless,  we  can  show  the  basis  function  Wa  approaches  the 
basis  function  for  the  TPS  as  a  —  Q.      Since  the  Green's  function 
for  the  thin  plate  (the  biharmonic  equation)  is  a  constant  times 
r  log  r,  which  changes  sign  at  r  =  1,  it  is  not  obvious  that  this 
will  happen.   In  order  to  show  that  this  is  the  case,  an  expan- 
sion of  KQ,  KQ(ar)  =  -(log(ar)  +  Y)(l  +  (ar)2/4  +  0((ar)2)  for 
small  values  of  <*  (see  Abramowitz  and  Stegun  (1964),  p.  375)  is 
substituted  into  the  expression  (2).   Here  Y  is  Euler's  constant. 
Upon  integration  of  the  first  few  terms,  one  obtains 

W  (r)  =  -(2Tr)"1r2log  r  +  (r2/(8iT))(l  -  log(a/2)  -Y)  +  0(«2). 

Since  additive  r   terms  in  the  basis  do  not  change  the  TPS  (it  is 
invariant  with  respect  to  the  units  used  to  measure  distance),  it 
is  seen  that  as  a — -0,  the  basis  functions  do  approach  those  of 
the  TPS. 


The  TPST  is  also  invariant  under  change  of  the  units  used  to 
measure  distance,  provided  the  tension  parameter  is  also  changed 
to  the  new  units.   If  r  is  replaced  with  fir,    and  <*  with  a/p   ,    a 
change  of  the  variables  of  integration  in  (2)  to  obtain  the  same 

integral  form  reveals  that  the  basis  function  is  multiplied  by  a 

2 
constant,  fi    . 

Before  considering  the  variational  formulation  of  the  TPST, 

the  properties  of  the  approximation  as  r — -00  will  be  investigated. 

Consideration  of  Wa  (r)  for  large  r  (r  >  R)  shows  that 

r   t 

Wa(r)=Wa(R)  -  RW^(R)log  R  +  RW^(R)log  r  +  (2 v) 'lft'1f sKQ (as) dsc 

o   o 

The  asymptotic  approximation   for    large  argument    (see  Abramowitz 

and  Stegun    (1964),    p.    378)    is 

K0(as)  ~  (V(2as))~1/2exp(-as)  (1  +  Ofs""1))  .   Substitution  of  this 

into  the  integral  above  shows  it  is  asymptotic  to 

(  V(2a3))-1/2(R1/2exp(-«R)log  r  +  CR  +  0(exp(-ar)), 
where  CR  is  a  constant  depending  on  R.   It  is  then  seen  that 

Wa(r)  =  DR  +  ER  log  r  +  0 (exp (-ar ) ) , 
where  DR  and  ER  are  constants  depending  on  R.   Thus 
Wa(r)  =  O(log  r).   The  interpolant  (3)  is 


N 


F(x,y)  =  TAkwa(rk>  +  a'   witn    / 
k^L  k=l 


Ak    =    0. 


Using    r^   =    (x  -    xk)2   +    (y  -    yk)2  =  r2  -    2P|crcos9k  +   Pk   = 
r2(l    -    2(pk/r)cos9k    +    Pk/r2)    for    large    r,    one    obtains 


F(x,y)  =  V  Ak(DR+(ER/2)log  [r2 (1-2 (  Pk/r)cosek+  Pj*/r2)  ] 
k  =  l 


+   a   +  0(exp(-ar) ) 


N 


a  ,-i 


=    (ER/2)  \   Ak  [log    r^+log(l-2(  Pk/r)cos0k+ P^/r2)] 


+   a   +  O(exp(-or) ) 
=    (ER/2)  V  Ak[-2(pk/r)cos9k+P^/r2+0(r"2)] 


N 

I 


+   a   +   0(exp(-ar) ) 


=   a   +    OU"1)     . 


Using  the  same  ideas,  it  is  easy  to  show  that  the  first  partial 
derivatives  of  f  are  0(r   )  as  r— °°. 

Consideration  of  the  variational  form  of  the  TPST  gives  the 
minimization  property  for  these  approximations.   The  functional 
is  the  expected  generalization  of  that  for  univariate  spline 
under  tension,  given  by 


.2     k  2 


rr      /3w\*        o/l_W_  \        l/3w\x      2  /  3w  \        .       2  /  3W  \ 


2   2 
_v 


The  derivation  of  the  above  yields  an  additional  term  which 

depends  only  on  the  data.   Proof  that  this  functional  has  finite 

_  j 

value  relies  on  the  fact  that  the  first  derivatives  are  0(r   ), 

going  to  zero  sufficiently  fast  for  the  integrals  to  exist. 

The  above  approach  yields  a  method  with  polynomial  precison 
for  constants  only.  If  the  data  all  lie  on  a  plane  not  parallel 
with  the  xy-plane,  the  resulting  surface  will  not  be  a  plane. 


This  arises  from  the  application  of  tension  in  the  horizontal 
direction.   It  is  trivial  to  obtain  linear  precision  by  inclusion 
of  the  two  linear  terms,  as  indicated  in  (1).   The  physical 
analogue  of  this  is  not  immediately  apparent,  but  the  resulting 
surfaces  are  more  consistent  with  the  surfaces  obtained  using 
TPS,  and  perhaps  more  satisfactory  overall.   Of  course,  the  above 
functional  cannot  be  finite  if  the  linear  terms  are  included  in 
the  approximation.   We  give  examples  of  approximations  computed 
both  ways  in  the  next  section. 

4.   Numerical  Examples 

Three  examples  will  be  given  which  exhibit  a  variety  of 
behavior,  both  good  and  bad.   All  of  the  examples  are  defined  on 
the  unit  square,  so  the  tension  values  in  the  three  cases  have 
the  same  units.   With  the  exception  of  Figure  5,  all  approxima- 
tions were  computed  with  the  inclusion  of  linear  terms.   Figure  5 
was  computed  with  only  the  constant  term  for  purposes  of  showing 
a  comparison. 

The  first  example  is  a  series  of  surfaces  which  all  interp- 
olate the  same  data,  but  with  increasing  tension.   The  data  is 
from  Nielson  and  Franke  (1984,  Table  2).   There  are  23  points. 
No  tension  results  in  overshoot  in  the  vicinity  of  the  sharp 
gradient.   Figure  2  shows  that  as  the  tension  is  increased,  the 
overshoot  becomes  less,  and  finally  disappears.   In  this  example 
the  surface  is  generally  well  behaved  as  tension  is  increased. 

A  second  example  is  also  from  the  paper  by  Nielson  and 
Franke  (Table  3)  and  involves  only  eight  points.   Figure  3  shows 
the  emergence  of  the  sharp  "points"  of, the  basis  functions  at  the 


10 


data  points  as  tension  increases,  while  away  from  the  data  the 
function  remains  very  smooth,  in  contrast  to  the  piecewise  linear 
function  over  the  tr iangulation  approached  by  the  surface  in  the 
reference.   The  behavior  appears  to  be  somewhat  like  one  would 
imagine  a  rubber  sheet  would  behave:   too  thick  to  be  a  membrane, 
but  supporting  little  displacement  except  near  the  data. 

The  third  example  shows  that  while  the  use  of  tension  may 
help  to  control  overshoot,  it  may  lead  to  undesirable  behavior 
away  from  the  sharp  gradient.  The  data  is  taken  from  Franke  and 
Nielson  (1983,  Table  1),  has  33  data  points,  and  is  an  extreme 
example  since  the  parent  surface  is  discontinuous.   Figure  4 
shows  that  as  the  tension  is  increased  the  approximation  tends  to 
look  a  bit  like  a  plane  with  local  sharp  peaks  and  dips  to 
achieve  interpolation.   Thus,  while  moderate  values  of  tension 
result  in  improvement  of  the  overshoot,  as  tension  increases  the 
effect  of  this  localization  degrades  the  surface.   For  compar- 
ison, Figure  5  was  computed  with  the  approximation  including  only 
the  constant  term.   In  part  (a),  where  the  tension  has  a  moderate 
value,  little  difference  is  seen.   In  part  (b),  for  a  larger 
value  of  tension  the  character  of  the  surface  on  the  back  "flat" 
part  is  seen  to  be  very  different,  with  the  surface  being  pulled 
down  to  the  data  in  Figure  4,  while  being  pulled  up  to  the  data 
in  Figure  5. 

Other  surfaces  have  been  computed  including  only  the  con- 
stant term,  in  particular,  all  of  the  surfaces  illustrated  above. 
The  differences  between  them  are  small,  for  the  most  part,  the 
example  given  in  Figure  5  being  the  extreme  case  investigated. 
The  additional  examples  are  shown  in  the  appendix. 
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5.   Conclusions 

The  dispacement  of  a  thin  plate  with  tension  has  been  devel- 
oped and  applied  to  the  problem  of  scattered  data  interpolation. 
It  appears  the  scheme  may  give  the  user  a  useful  way  of  control- 
ling the  behavior  of  the  interpolating  surface.   In  one  variable 
some  work  has  been  accomplished  toward  determination  of  a  suit- 
able value  of  the  tension  parameter  (see  Lynch  (1982)).   In  the 
present  situation,  there  is  no  known  way  to  automatically  choose 
the  tension  parameter  to  achieve  the  proper  amount  of  control. 
However,  this  is  not  in  marked  contrast  to  the  whole  problem  of 
scattered  data  interpolation,  where  one  should  proceed  cautiously 
until  it  is  determined  how  one  or  more  methods  perform  for  data 
sets  with  the  given  configuration  and  features.   Other  issues, 
such  as  scaling  of  the  data  may  alter  the  approximation  in  a 
significant  way  (e.g.,  see  Breaker  (1983)),  and  may  be  important 
here  as  well. 

The  work  of  Duchon  resulted  in  the  characterization  of  the 
TPS  in  terms  of  minimization  of  the  thin  plate  functional  in  a 
certain  reproducing  kernel  Hilbert  space.   It  is  almost  certain 
that  the  TPST  also  minimizes  the  energy  functional  noted  above  in 
a  similar  space. 
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Figure  3 
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Appendix 

This  appendix  is  given  to  illustrate  more  fully  the  behavior 
of  the  TPST  as  well  as  to  document  the  differences  between  the 
two  variants  of  it.   The  main  content  is  a  set  of  full  page  plots 
corresponding  to  a  number  of  different  data  sets  interpolated 
with  various  tension  values,  using  the  linear  terms,  on  the  left 
(even  numbered  pages,  Figure  n.a),  and  the  same  data  interpolated 
using  only  the  constant  term,  on  the  right  (odd  numbered  pages, 
Figure  n.b).   This  enables  one  to  make  quick  comparisons  between 
the  two  surfaces.   Table  1  is  a  table  of  contents  for  the  plots, 
identifying  the  Figure  number  in  which  the  surface  appears.   For 
completeness  I  have  included  tables  of  the  data,  in  Tables  2-6. 

It  will  be  seen  that  there  are  only  small  differences  in  the 
corresponding  surfaces  for  small  values  of  tension  (<*  =  10).   In 
some  cases  the  surfaces  are  very  similar  for  large  values  of 
tension  (  <*  =  50)  as  well,  although  in  other  cases  the  surfaces 
are  significantly  different,  as  in  the  surfaces  shown  in  Figures 
4  and  5.   These  same  surfaces  are  shown  in  Figures  15a  and  17a  in 
this  appendix,  and  it  is  seen  that  the  more  pleasing  of  the  two 
is  the  surface  including  only  the  constant  term.   On  the  other 
hand,  perusal  of  the  surfaces  shown  in  Figures  13  and  14  shows 
the  more  pleasing  surface  to  be  that  including  the  linear  terms. 
Figures  13a  and  14a  are  the  same  surface  as  shown  in  Figure  2d 
and  2f.   The  other  examples  may  show  one  or  the  other  to  be  more 
pleasing;   to  a  certain  extent  it  is  a  personal  choice. 

It  is  apparent  that  neither  of  the  schemes  will  be  the 
better  in  all  cases.   There  are  good  reasons  for  wanting  such  a 
method  to  have  linear  precision,  and  I  have  favored  the  scheme 
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incorporating  the  linear  terms.   Nonetheless,  in  certain  applica- 
tions it  may  well  be  advantageous  to  include  only  the  constant 
term,  and  this  is  easy  to  justify.   Thus,  the  user  can  decide 
based  on  individual  mathematical  preference,  or  on  which  scheme 
gives  the  more  pleasing  surface. 
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0.0000 
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9125 
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0.0000 
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7125 

0 

7625 

0.0000 
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2250 
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0.0000 

0 

0500 

0 
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0.0000 

0 

4750 

0 

6375 

0.0000 

0 

8500 

0 

4375 

0.0000 

0 

7000 

0 

3125 

0.0000 

0 

2750 

0 

4250 

0.0000 

0 

4500 

0 

2875 

0.0000 

0 

8125 

0 

1875 

0.0000 

0 

4500 

-0 

0375 

0.0000 

0. 

0500 

-0 

0500 

0.0000 
Table  1: 

0 

4500 

1 

0375 

0 

0000 

1 

0875 

0 

5500 

0 

0000 

0 

5375 

0 

8000 

0 

0000 

0 

0375 

0 

7500 

0 

0000 

0 

1875 

0 

5750 

0 

0000 

0 

7125 

0 

5500 

0 

0000 

1 

0000 

0 

2625 

0 

0000 

0 

5000 

0 

4625 

0 

0000 

0 

1875 

0 

2625 

0 

2000 

0 

5875 

0 

1250 

0 

0000 

1 

0500 

-0 

0613 

0 

0000 

0 

1000 

0 

1125 

0 

0000 

CARDINAL  FUNCTION 


0.0000 
0.0000 
0.3000 
0.8000 


0.0000 
1.0000 
0.7000 
0.8200 


0.4000 
0.2000 
0.4000 
0.1000 


1.0000 
0.2000 
0.8000 
1.0000 


0.0000 
0.3000 
0.2000 
1.0000 


Table  2:   EIGHT  POINT 


0.3000 
0.4000 
0.4000 
0.0000 
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0.0000 

0.0000 

0.3000 

0.0000 

0.6000 

0.3600 

0.0000 

1.0000 

0.3000 

0.2400 

0.0400 

0.2910 

0.1600 

0.3200 

0.3000 

0.2000 

0.6000 

0.3300 

0.2000 

0.8000 

0.3000 

0.2000 

1.0000 

0.2400 

0.4000 

0.4000 

0.3000 

0.3600 

0.5600 

0.3000 

0.3200 

0.7600 

0.0750 

0.4000 

1.0000 

0.1500 

0.4800 

0.4800 

0.3000 

0.5200 

0.5200 

0.0750 

0.5200 

0.6400 

0.0750 

0.6000 

0.0800 

0.0210 

0.6000 

0.4000 

0.2250 

0.6000 

0.8000 

0.0750 

0.6800 

0.5600 

0.0750 

0.8000 

0.2000 

0.1500 

1.0000 

0.0000 

0.1200 

0.9600 

0.6000 

0.0600 

1.0000 

1.0000 

0.0000 

Tabl 

e  3:   ROCKY 

MTN  JRNL  FNCTN 

X 


X 


0.0400 

0.0400 

0.5000 

0.4000 

0.0400 

0.5000 

0.6000 

0.0400 

0.5000 

0.8000 

0.0000 

0.5000 

0.9600 

0.0800 

0.5000 

0.0000 

0.2800 

0.5000 

0.2000 

0.2700 

0.5000 

0.4000 

0.2000 

0.5000 

0.8800 

0.2000 

0.5000 

0.1600 

0.4000 

0.5000 

0.2800 

0.3800 

0.5000 

0.4800 

0.3600 

0.5000 

0.6800 

0.3200 

0.5000 

.  0.8000 

0.3700 

0.5000 

1.0000 

0.3600 

0.5000 

0.4000 

0.4400 

0.3267 

0.6400 

0.4400 

0.1960 

0.8400 

0.4200 

0.0934 

1.0000 

0.4400 

0.0000 

0.7200 

0.5600 

0.0941 

0.9200 

0.6000 

0.0222 

0.4800 

0.6400 

0.1170 

0.9600 

0.7600 

0.0040 

0.6400 

0.7600 

0.0360 

0.4000 

0.8800 

0.0150 

0.8000 

0.8800 

0.0050 

0.6000 

1.0000 

0.0000 

1.0000 

1.0000 

0.0000 

0.0800 

0.4800 

0.4911 

0.0000 

0.6000 

0.4444 

0.0400 

0.8000 

0.2778 

0.0000 

1.0000 

0.0000 

0.2800 

1.0000 

0.0000 
Table  4: 

FAULT  LINE  SIX 
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X 


X 


0.0400 
0.6000 
0.9600 
0.2000 
0.8800 
0.2800 
0.6800 
1.0000 
0.6400 
1.0000 
0.9200 
0.9600 
0.4000 
0.6000 
0.0800 
0.0400 
0.2800 
0.0540 
0.0175 
•0.0510 
•0.0487 
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0.1090 
0.0935 
0.1453 
0.0627 
0.2392 
0.2767 
0.2267 
0.2426 
0.3858 
0.3179 
0.3777 
0.3813 
0.2804 
0.4278 
0.4664 
0.4092 
0.4812 
0.4027 
0.5730 
0.5014 
0.6107 
0.5381 
0.5026 
0.6428 
0.6704 
0.6334 
0.6896 
0.6838 
0.7635 
0.8259 
0.8087 
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0.8171 
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0.8600 
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9471 


0 
0 
0 


1 
0 


0.0400 
0.0400 
0.0800 
0.2700 
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0.3800 
0.3200 
0.3600 
0.4400 
0.4400 
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0.8800 
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1 
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0 
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0.2212 
0.0005 
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0.2490 
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0.0010 
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0.2478 
0.0558 
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0.0008 
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0. 
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0. 
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0. 
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0. 
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0. 
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0 
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0. 
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0. 
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0. 
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0. 
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1. 
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1. 

0000 

0. 
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0. 

0000 

0. 
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0. 
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0 
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0. 

0000 

0. 
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0. 
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0 
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0. 
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0 
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0 

0395 

0. 
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0. 
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0 
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0. 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 

0768 

0 

4171 

0 

4996 

0 

0959 

0 
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0 

2646 

0 

0293 

0 

5000 

0 

2089 

0 

2669 

0 
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0 

1715 

1 

0000 

0 

0000 

0 
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0 

9046 

0 

0122 

0 

3663 

0 

0397 

0 

5000 

0 

3832 

0 

2390 

0 
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0 
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0 
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6445 

0 
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0 
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0 
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0 
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-0 

0285 
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5000 
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4200 

0 
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5000 
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0 
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0 
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0 

4793 

0 
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0 
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0 
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0239 

0 
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-0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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-0 
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0 
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0 
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0 
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0 
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-0 
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1 
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0 
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-0 
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1 

0415 

0 
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-0 

0012 

1 

0000 

0 
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0 
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Data  Set  Table 
Cardinal  1 
8  Point  2 
Rocky  Mtn  J  F  3 
Fault  Line  6  4 
Fault  Line  7     5 


Tension  =  <* 

(=TN) 

10 

30 

50 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Table  6:   Figure  number  in  which  surface  plot  appears 
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25  PT   CARDINAL  FUNCTION 
THIN  PLATE  WITH  TENSION 

Figure    6. a 


61984 
TN    = 


1437 
10.0 
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25  PT   CARDINAL  FUNCTION 
TPS  WITH  TENSION,  CONST 

Figure    6.b 


121084 
TN    = 


1655 

10.0 


25 


25  PT   CRRDINRL  FUNCTION 
THIN  PLATE  WITH  TENSION 


121084 
TN    = 


1604 
30.0 


Figure    7- a 
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25  PT   CARDINAL  FUNCTION 
TPS  WITH  TENSION,  CONST 


121  184 
TN    = 


1  158 
30.  0 


Figure    7.b 
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25  PT   CRRDINRL  FUNCTION 
THIN  PLRTE  WITH  TENSION 


Figure    8. a 


61984 
TN    = 


1437 
50.0 
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25  PT   CARDINAL  FUNCTION 
TPS  NITH  TENSION,  CONST 
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Figure    8.b 
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8  PT   EIGHT  POINT 
THIN  PLATE  WITH  TENSION 
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Figure    9. a 
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8  PT    8  POINT 
TPS  WITH  TENSION 
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Figure    9«b 
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8  PT   8  POINT 
THIN  PLATE  WITH  TENSION 
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Figure    10. a 
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8  PT   8  POINT 
TPS  WITH  TENSION, 
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Figure    10. b 
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8  PT   EIGHT  POINT 
THIN  PLRTE  WITH  TENSION 
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Figure    11. a 
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8  PT 
TPS  NITH 
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TENSION, 
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Figure    11. b 
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23  PT   ROCKY  MTN  JRNL  FNCTN 
THIN  PLRTE  NITH  TENSION 
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Figure    12. a 
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23  PT   ROCKY  MTN 
TPS  NITH  TENSION, 


JRNL  FNCTN 
CONST 
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Figure   12. b 
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23  PT   ROCKY  MTN  JRNL  FNCTN 
THIN  PLRTE  WITH  TENSION 


62184 
TN    - 


827 

30.0 


Figure    13. a 
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23  PT   ROCKY  MTN  JRNL  FNCTN 
TPS  WITH  TENSION,  CONST 
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Figure    13. t> 
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23  PT   ROCKY  MTN  JRNL  FNCTN 
THIN  PLATE  WITH  TENSION 
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Figure   l4.a 
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23  PT   ROCKY  MTN 
TPS  WITH  TENSION, 


JRNL  FNCTN 
CONST 

Figure    l4.b 
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33  PT   FAULT  LINE  SIX 
THIN  PLATE  WITH  TENSION 
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Figure   15. a 
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33  PT   FAULT  LINE  SIX 
TPS  WITH  TENSION,  CONST 
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Figure   15. b 
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33  PT   FRULT  LINE  SIX 
THIN  PLATE  NITH  TENSION 
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Figure    l6.a 
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33  PT   FRULT  LINE  SIX 
TPS  WITH  TENSION,  CONST 
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Figure    l6.b 
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33  PT   FAULT  LINE  SIX 
THIN  PLATE  NITH  TENSION 
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Figure   17. a 
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33  PT   FAULT  LINE  SIX 
TPS  WITH  TENSION,  CONST 
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Figure   17. b 
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Figure    l8.a 
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130  PT   FAULT  LINE  SEVEN 
TPS  WITH  TENSION,  CONST 
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Figure    l8.b 
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130  PT   FRULT  LINE  SEVEN 
THIN  PLRTE  NITH  TENSION 
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TN  =       30.0 


Figure  19. a 
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130  PT   FAULT  LINE  SEVEN 
TPS  NITH  TENSION,  CONST 
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Figure    19. b 
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130  PT   FAULT  LINE  SEVEN 
THIN  PLATE  WITH  TENSION 
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Figure    20. a 
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.30  PT   FAULT  LINE  SEVEN 
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Figure    20. b 
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